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Abstract 

We introduce multilinear operators, that generalize Hirota's bilinear D operator, based 
on the principle of gauge invariance of the r functions. We show that these operators can 
be constructed systematically using the bilinear D^s as building blocks. We concentrate 
in particular on the trilinear case and study the possible integrability of equations 
with one dependent variable. The 5th order equation of the Lax-hierarchy as well 
as Satsuma's lowest-order gauge invariant equation are shown to have simple trilinear 
expressions. The formalism can be extended to an arbitrary degree of multilinearity. 
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1. Introduction 



The Hirota bilinear operators were introduced as an antisymmetric extension to the 
usual derivative [1], because of their usefulness for the computation of multisoliton 
solution of nonlinear evolution equations. The bilinear operator Dy; = — d^^ acts on 
a pair of functions (the 'dot product') antisymmetrically: 

Dxf-g = {dx, - dx,)f{xi)g{x2)l^^^^^^ = f'g - g'f. (1) 

The Hirota bilinear formalism has been instrumental in the derivation of the multisoliton 
solutions of (integrable) nonlinear equations. A prerequisite to its application is a 
dependent variable transformation that converts the nonlinear equation into a quadratic 
"prepotential" form. This is best understood in a specific example, so let us consider 
the paradigmatic case of the KdV equation. Starting from 

Uxxx + Quux + Wt = 0, (2) 

we introduce the transformation u = 2dl logF and obtain (after one integration): 

FxxxxF — 4:FxxxFx + SF^x + FxtF — FxFt = 0. (3) 

This last equation can be written in a particularly condensed form using the Hirota D 
operator: 

{Dt + DxDt)F-F = 0. (4) 

The power of the bilinear formalism lies in the fact that for multi-soliton solutions the 
F's are simple polynomials of exponentials [2] . Thus the construction of soliton solutions 
becomes an algebraic problem. This approach has made possible the investigation of 
large classes of bilinear equations and the classification of integrable cases [3]. (The 
integrability of these equations has been confirmed by singularity analysis [4]). 

2. Gauge-invariant bilinear operators 

An important observation (that has motivated the present line of research) is the re- 
lation of the "physical" variable u = 2dx log F (for the KdV equation) to the Hirota's 
function F: the gauge transformation F — > e^^+'^*F leaves u invariant. It turns out 
that this is a general property of bilinear equations. In fact, one can deGne the Hirota's 
bilinear equations through the requirement of gauge invariance. We will now prove this 
statement. 

Let us introduce a general bilinear expression and ask that it be invariant under 
the gauge transformation F ^ e^F with 77 = px -\- ut: 



N N 

Y,Ck{d^ey){d''-^e^g) = e^^ J] Cfe(aV)(5^-'^^). 

k=0 k=0 



(5) 



Expanding the left hand side and equating the coefficients of {d''f) {d^ ^g) we obtain, 
for all n, m: 

Ef!) ("""Mp^-^c.^c.^..-™ (6) 



fc=0 



where f ^ | is the binomial coefficient. From the structure of the left hand side, we 

have the inequalities 0<n<k<N — m. Taking m = N — n we find that k = n and 

(6) is identically satisfied. For all other values of m (m =0, 1, . . . , N — n — 1) we must 
find those Cfc's that satisfy 

In the special case m = N — n — 1 (7) reduces to: 

{N - n)cn + {n+ l)cn+i = 0, 

the solution of which is 

Cn = (-ir(^)co. (8) 

This solution does, in fact, satisfy (7) for all values of m: using (8) we find that the Ihs of 

(7) is just the expression of (1 — l)-^~"^~"^Co and thus equal to zero. So the only gauge in- 
variant bilinear differential operator is (up to a scaling of cq) X^^o(~^)'^ ^ ) did2~^ = 

(di — 82)^ 1 i.e. the Hirota operator (The use of the indices (1,2) may appear su- 
perfluous at this stage, since in the bilinear case there are only two variables on which D 
operators act. However the notation will be particularly useful in the higher multilinear 
cases) . 

3. Gauge-invariant trilinear operators 

The Hirota bilinear operators can be obtained on the sole requirement (admittedly a 
strong one) of gauge invariance. In this paper our objective is the extension of the 
bilinear formalism and the introduction of multilinear operators. Very few results ex- 
ist in this direction. Satsuma and collaborators have introduced a particular class of 
trilinear equations, that can be written as a single (3 x 3) determinant [5]. A full hier- 
archy of equations was obtained and the richness of the solutions presented are a strong 
indication of their integrability. However, Satsuma's approach offers no clue on what 
a trilinear operator should be. Here we will use the same principle used to find the 
bilinear L>'s: gauge invariance. 

For trilinear expressions the invariance condition writes: 

ckimid'ey)id'e^g)id^e^h) = e^^ J2 cum{d'' f){d'g){d^h). (9) 

k+l+m=N k+l+Tn=N 



In analogy to the bilinear case we find: 



k+l+m=N ^ ^ ^ / \ ^ / 

For K + \ + jj, = N equation (10) is identically satisfied and we are left with 

S "'^'-f^) (a) ^^'^^ « + A + /x<N). (11) 

Consider the k; + A + // = A'' — 1 equations. They write: 

{jj. + l)c^AAt+i + (/« + 1)ck+iam + (A + l)c«A+iM = 0- (12) 

We have ^^^^'^^ such equations for the (i^[±lK^^±^^ ^'s. Thus these equations would 
determine the c's up to A?" + 1 free coefficients, provided the rank of the system is 
maximal. This is indeed the case. In fact, the equations (k. A, = 0) contain for the 
first time c«ai (which can be expressed in terms of the (AT + 1) c^ao)- We can then 
compute successively the higher-// terms up to the last equation (k = 0, A = 0, // = N—1) 
from which we can solve for cqoat-i- Thus all c^\^j_ for > can be expressed in terms 
of the c^Ao- So, given that the rank is maximal, we can choose any basis for the c's. A 
most convenient basis are the following N -\- 1 operators: {di — d2)^{di — Ss)"^""^ for 
n = 0, ...,A^. 

Thus the basic building blocks for the trilinear operators are again the Hirota 
bilinear D's, we must just specify the indices in this case. We thus have D12 = d^^ —d^^, 
D23 = —dx3, -D31 = —dx^, but, of course the three are not linearly independent: 
-D12 + -D23 + -D31 = 0. Their action on a 'triple dot product' is analogous to the bilinear 
case: 

D^2f-9-h = {dx, - ^x,)fix^)g{x2)h{xs)\^^^^^^^^^^ = {fg - fg')h. (13) 

The choice of a particular pair of D's as the basic trilinear operators breaks the sym- 
metry between the three coordinates x^'s. It is possible to restore this symmetry by 
introducing a different basis for the trilinear operators, T and T*: 

T = d^+ jd2 + fds , T* = ai + fd2 + jds, (14) 

where j is the cubic root of unity, j = e^*'^/^. (Note that the star in T* indicates 
complex conjugation for the coeflBcients in T but not for the independent variables). 
The price we have to pay for restoring this symmetry is that the operators are now 
more complicated. Note that T^T*^F-F-F = unless n - m = (mod 3), which is 
the equivalent to the bilinear property D^F-F = unless n = (mod 2). 

The generalization to higher multilinear equations is straightforward. One can 
introduce the set of n{n — l)/2 operators Dij acting on n-tuple dot-products -Dy/i-/2- 
. . . ■ fn- (Of course only n — 1 of the D^'s are independent, a convenient basis being 



the Dij, j = 2, . . .n). As in the trihnear case, one can also construct "symmetric" 
operators: 

n-l 

Mm^Yl ^'mdk (15) 
k=l 

where the Zm^s are the (n — 1) n-th roots of unity other than one. 

4. Examples of multilinear equations 

Multihnear operators are not just a trivial extension of the Hirota bilinear formalism. 
They are necessary for the description of nonlinear evolution equations that cannot 
be cast in a bilinear form and such equations do exist. An interesting example is the 
fifth-order equation of the Lax hierarchy [6]: 

Uxxxxx + lOuUxxx + ^OUxUxx + ^Ou^Ux + Uf = 0. (16) 

While this Lax-5 equation does not possess a simple bilinear expression like KdV itself, 
it has a trilinear form (with u = 2d1 log F) : 

(7T| + 20T^Tf + 21TxTt)F-F-F = 0. (17) 

In the previous section we referred to Satsuma's trilinear equations [5] . The lowest-order 
one 

F F F 

(18) 



(equivalent, through F = e"', to the Monge- Ampere equation w^y — WxxWyy = 0) is 
gauge-invariant and can be written as: 

{TxT*TyT* - TlTf)F-F-F = 0. (19) 

(The higher Satsuma equations are given by sets of equations with "dummy" indepen- 
dent variables and it is not clear how to implement the gauge-invariance requirement in 
such a situation). 
The equation: 

p'^p — fPfP fP —AFFF +2FF F A- AF F F — 2F F"^ 
^ ^ xxxxy ^ ^ xxxx ^ xxxy T -^-t -t xx^ xxy \ x^ y^ xxx y^ xx 

-AFxFxxFxy + ^FlFxxy + ^{F'' Fxxt - FFxxFt - 2FFxFxt - 2F^Ft) = 0, (20) 

obtained as a reduction of a self-dual Yang- Mills equations [7] , can also be written as a 
genuinely trilinear equation: 

(T^T* + 8T^T*Ty - S6T^Tt)FFF = 0. (21) 
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Further examples can be presented and we can, of course, construct also higher 
multilinear equations (quadri-, penta-, etc.) equations. Instead of dealing with specific 



cases let us present here some general considerations. Let us start with a nonlinear 
(in u) equation, of order k + 1 having the form ut + dxP{u, Ux, ■ ■ ■ , Ukx) = 0. Several 
well known integrable equations belong to this class. It is simpler to work with the 
time independent part Ux, . . • , Ukx) = which only involves derivatives up to order 
k. We consider the leading part of P which we assume to be weight-homogeneous in 
u and dx-i with u having the same weight as d^. Then we can transform this leading 
part to a multilinear expression through the transformation u = ad1{\ogF) and obtain 
generically an /c + 2-multilinear equation. The scaling factor a can then be chosen 
so as to make the lowest-derivative terms vanish, that is those that appear under the 
combination _ h^p" p'kp'^_ (This is possible because these two terms have 

a common factor, polynomial in a). Then an F"^ term can be factored out and the 
resulting multilinear equation is at most /c-linear. 

At order five {k = 4) we have three integrable equations. We should expect, in 
principle, these equations to have quadrilinear forms. Some unexpected cancellations, 
however, do occur. Thus the Sawada-Kotera [8] equation has a bilinear expression, the 
Lax-5 has the trilinear form we gave in (17), but the Kaup-Kuperschmidt [9] is quadrilin- 
ear. At order seven [k = 6) three integrable equations were known to exist. The higher 
Sawada-Kotera has a trilinear expression (see next section) while the 7th order equation 
in the Lax hierarchy is a pentalinear one. Again, for the higher Kaup-Kuperschmidt 
equation no extra simplification is possible and this equation is hexalinear. For all these 
equations the time derivative can be incorporated in the multilinear equations in a very 
simple way without altering the degree of multilinearity. 

Thus the multilinear extension to Hirota bilinear approach has a wide range of 
applicability (in particular if we allow for multicomponent equations, introducing more 
than one dependent functions, in analogy to the bilinear case). 

5. Singularity analysis of trilinear equations 

All the above equations have as a common characteristic their integrability. The study 
of integrability is, in fact, the motivation behind the multilinear approach. The sys- 
tematic classification of bilinear equations we presented in [3,4] was based on the study 
of multisoliton solutions and of the Painleve property. We intend to come back to the 
investigation of soliton solutions for our multilinear equations in some future work. In 
the following paragraph we will limit ourselves to the singularity analysis of trilinear 
equations involving only one dependent variable, i.e. unicomponent equations. (Let 
us recall here that in the bilinear case the study of these simplest equations led us to 
conclusive results). 

In order to perform the Painleve analysis, we shall study the leading (highest- 
order) part of the equations with just one independent variable. This is sufficient for 
the computation of dominant singularities and resonances, although for the check of 
resonance compatibility we would need the full equation, which remains unspecified at 
this stage. (For example, for the analysis of an equation like (17), we would consider 



only the 7^6 ji3ji*3 terms and not T^Ti^ since the influence of the latter would appear 
only at the resonance condition). 

Since the dependent function m in a nonlinear equation is related to the multilinear 
F through u = log F it is clear that a zero in F induces a pole- like behaviour in 
u. Let us show how one performs the singularity analysis for trilinear equations in a 
concrete example: T^T*^ f ■ f ■ f — 0. Putting / ~ x"^ (for the dominant part) we find 
n = 0, 1 as the only possible behaviours. The first corresponds to a nonsingular Taylor- 
like expansion, which is always possible. The second behaviour / ~ x corresponds to 
a simple zero that would give a (double) pole in u. The resonances in this case are 
obtained if we substitute / = x{l + (px"^) and collect terms linear in (p. The result 
is r = —1,0, 1,6 and (after a check that no incompatibilities arise at any resonance) 
we conclude that this equation passes the Painleve test. This is what one would have 
expected, had we looked at the nonlinear form of the equation, u^x + = 0, which is 
just the time-independent part of KdV integrated once. 

We shall not present the details of the singularity analysis of all the equations 
that we have studied. The results are summarized below. The notation we are using 
is the following: E{n,m), (which is identical to E{m,n)), represents the expression 
j'mrp*np_p_p Note that for a given N = n + m there may exist several pairs of (n, m) 
such that E{n,m) is not identically zero, namely those for which n = m (mod 3). The 
leading part of the general equation at order N is then given by a linear combination 
of all the non- vanishing E{n,mys. In each case we give below the precise combinations 
that lead to equations with the Painleve property. The nonlinear forms of the equations 
are obtained by the standard substitution F = followed by tt = 2g". 
N = 2: E{1,1) 

In this case we can write the result also in bilinear form: 

E{1, 1) oc F{D'^F-F) = 2F{F"F - F'^) = 2e^3g" = e^9u. 
N = S: E(S,0) 

E{3, 0) oc F'"F2 - 3F"F'F + 2F'^ = e^^g'" oc e^^u'. 
iV = 4: F(2,2) 

Here also we can write the result in bilinear form: 

E{2, 2) oc FiD'^F-F) = 2e^9{g"" + Qg"'^) = e^^{u" + Sw^). 

This, of course, is just the leading part of the KdV equation in potential form. 
iV = 5: F(4,l) 

F(4,1)/F3 oc'u'" + 6m'm. 

Note that this is the derivative of the expression obtained at A?" = 4. 
N = ^ : AF(6,0) + |Ui?(3,3). This is the first case where we have two possible n, m 
pairs. The A, \x combinations that pass the Painleve test are the following 

a) (7F(6, 0) + 20F(3, 3))/^^ oc v!'" + X^vil'u + 5w'2 + \^u^ , 

This is the leading part of the 5th order equation in the Lax hierarchy of KdV, 
eq.(16), integrated once. 

b) (-2F(6, 0) + 20F(3, 3))/^^ oc v!'" + Vbv!'u + 15^^. 



This is the leading part of the Sawada-Kotera equation, integrated once, 
c) (E(6, 0) - E{3, 3))/F3 oc uu" - u''^ + . 

The trihnear form of this case is F"" {F" F -F''^)-F"''^F + 2F'" F" F -F"^ , which 
can be cast in determinantal form and is a 1-dimensional member of the Satsuma 
family: 
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N = 7: E{5,2) 

E{5, 2)/F^ oc ^(5) + 15u"'u + 15u"u' + ^hu'u^. 

This is the Sawada-Kotera equation, i.e. the derivative of the expression obtained 
in AT = 6b. 

N = 8 : XE{7,1) + IJ,E{4, 4) and we take u = 6g" instead of tt = 2g" used before 

a) (4E(7, 1) + 5E(4, 4))/F3 oc u^^^ + 6u""u + 10u"'u' + 5u"'^ + lOw'^w^ + 10u''^u + | 
This would correspond to a 7*'* order equation which, we believe, leads to a new 
integrable case. 

b) {4:E{7, 1) + 14E(4, 4))/F3 oc + ^^''''u + lu"'v! + lu"'^ + \Au"u^ + Ivl'^u + 
This is the higher Sawada-Kotera equation we referred to in the previous section. 

c) (£;(7, 1) - £;(4, 4))/F3 oc v!"'u - ?>u"'u' + 2w"2 + ^^''u^ - Sw'^u + 

This new equation looks like an extension of Satsuma's equation given at N=6c 

above, but it cannot be written as a single determinant. 
At order = 9 there are no cases passing the Painleve test. 
N = 10: XE{8, 2) + fiE{5, 5) and we take u = 30g" 

{5E{8, 2) + 4E{5, 5))/F^ oc u^^^ + 2u^^\ + 4u^^^K' + 6u""u" + 5u"'^ + ^u""u'^ + 

Au"'u'u + 2u"^u + 2u"u'^ + ^u"u'^ + + 

This is also a new equation. 
Furthermore, no integrable candidates were found at orders N = 11, 12. The singularity 
analysis at these higher orders becomes progressively more difficult. (Already at AT = 12 
there exist three nonvanishing n,m combinations). It is, thus, not possible to extend 
our investigation to very high orders (as was done in our study of bilinear equations [4]), 
but we do believe that no further integrable candidates exist at higher orders. As in the 
bilinear finite (and rather small) number of unicomponent trilinear equations 

possess the Painleve property and can thus be integrable. 

6. Conclusion 

In the preceding paragraphs we have presented an extension of Hirota's bilinear for- 
malism that can encompass any degree of multilinearity. The main guide in our in- 
vestigation has been the requirement that the equations be gauge-invariant. Since our 
primary objective is the study of integrability, we have also presented a classification 
of one-component trilinear equations that pass the Painleve test. The crucial difference 
between the bilinear and the tri- (and multi-)linear case(s) is that now free parameters 
enter already at the leading part. This means that the Painleve analysis of the higher 



order unicomponent equations becomes increasingly diflBcult. Once the leading parts 
of these equations are fixed, one can study the lower-order terms that can be added 
without destroying the Painleve property. Starting from a complete classification of 
unicomponent equations one can build up multicomponent ones following the approach 
we presented in [3] for the bilinear case. Another interesting direction would be the 
computation of the multisoliton solutions of the trilinear equations. This would fur- 
nish another check for their possible integrability. Clearly, the domain of multilinear 
equations is still a terra incognita that deserves serious study. 
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